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Abstract 



We compare recent experimental results for one-particle-inclusive processes 
in deeply inelastic scattering at HERA with theoretical predictions in next- 
to-leading-order QCD perturbation theory, and study the factorization scale 
dependence of cross sections and charged multiplicities. In the future, for the 
HERA machine running at design luminosity, scaling violations of fragmenta- 
tion functions permit the measurement of the strong coupling constant. We 
estimate the size of the statistical error of a s that can be achieved, and study 
the theoretical error due to the various parton density parametrizations and 
due to the factorization scale dependence. 
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1 Introduction 



The DESY HI and ZEUS Collaborations have recently published results for 
charged-particle production in deeply inelastic electron-proton scattering at 
HERA In this paper, we calculate, in next-to- leading-order QCD per- 

turbation theory, the ^-distribution in the current fragmentation region, em- 
ploying the recent parametrization of charged-meson fragmentation functions 
from Refs. |||4[]- The comparison of theoretical and experimental results looks 
very encouraging and suggests that a determination of the strong coupling 
constant a s via scaling violations of fragmentation functions might be feasi- 
ble. We give an estimate of the statistical error Aaf at that can be expected 
for the HERA machine runnning at design luminosity, and study the errors 
Aaf DF and Aaf ale due to the uncertainty coming from various parton density 
parametrizations and choices of the factorization scale, respectively. Except for 
the dependence on the parton density parametrization, an a s -determination 
via scaling violations has the advantage that, being based on a renormalization 
group equation, it is model-independenQ A drawback is that the size of the 
effect, being proportional to the logarithm of the energy scale, is small. Com- 
pared with e + e~ annihilation, HERA offers the unique possibility to perform 
the required measurements at various energy scales at one machine without 
changing the centre-of-mass energy. 

The theoretical basis for the calculation of one-particle-inclusive cross sections 
is the factorization theorem of perturbative QCD (see, for example, Ref. |J 
and references therein). For the process IP —>■ I'h + X (I, I' are the incoming 
and outgoing leptons, P is the incoming proton, h is the observed hadron, and 
X denotes anything else in the hadronic final stat^j), the cross section a can 
be written in the form 

O- = / d£ f i/P (f , nf) / d^ D h/j (z, fijy) (Third (f , z, f%, Vd, f4) > (!) 



where fam and D h u denote parton densities and fragmentation functions, cr h J ard 
is the mass-factorized hard scattering cross section, £ and z are the momen- 
tum fractions of the incident parton i and of the observed hadron originating 
from the fragmenting parton j, and ///,//£> are the corresponding factorization 
scales. We have also indicated the renormalization scale fi r . The phenomeno- 
logical distribution functions / and D have to be taken from experiment. The 
hard scattering cross section <Thard has been calculated in Ref. 0. For the 



1 For example, in the case of (2+l)-jet production, the experimental jet rates have 
to be "corrected" back to the parton level by means of a fragmentation model, thus 
introducing an additional systematic uncertainty. 

2 We identify the four-momenta of particles with their genuine names. 



1 



present numerical study, we have used an implementation of a recent recal- 
culation 0. In leading order, the cross section is given by the process of the 
naive parton model with a fragmentation function attached to the outgoing 
quark, see Fig. |l|. 



Fig. 1. Feynman diagram corresponding to the leading-order contribution. 



In next-to-leading order, virtual corrections have to be added to the 7*g — > q 
diagram, and the QCD subprocesses 7*g — > qg and Y'g QQ have to be in- 
cluded. After factorization of collinear singularities in the parton cross section 
and renormalization of the phenomenological distribution functions, the re- 
sulting one-particle-inclusive cross section is infrared-finite. A comparison of 
theoretical and experimental ^-distributions is done in the next section. The 
scales fif, fiD an d fi r are in principle arbitrary. In the leading-order process un- 
der consideration, the only physical scale related to the hard scattering process 
is the photon virtuality Q = \/—q' 2 , with q the momentum of the exchanged 
photon. We will therefore identify the factorization and renormalization scales 
with Q, except for the case where we study the scale dependence explicitly. 

We will consider only processes with an exchanged virtual photon, and ne- 
glect the contributions from an exchanged Z boson. For the comparison with 
present-day experimental data, this is justified by the restricted range in Q. 
For simplicity, the analysis in Section |3] of the error of a possible measurement 
of the strong coupling constant via scaling violations of fragmentation func- 
tions is done in this approximation as well, even though Q may reach up to 
150 GeV. We assume this to be sufficient to achieve an estimate of the error, 
but it is clear that for an experimental analysis at large Q the contributions 
for an exchanged Z boson have to be included. 
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2 Comparison with Experimental Data 



In this section we compare the longitudinal momentum fraction distribution^] 
1 da 



p(x F ) 



(2) 



recently published by the HI || and ZEUS Collaborations with the the- 
oretical next-to-leading-order prediction. Here cr to t is the total cross section 
for the same cuts on the phase space of the outgoing lepton as for the one- 
particle- inclusive cross section a. The variable xp is defined to be 2hi,/W, 
where hi is the component of the observed hadron's momentum along the di- 
rection of the exchanged virtual photon in the hadronic centre-of-mass frame, 
(P + q) 2 is the total hadronic energy. The current direction is 
ion xp > 0. 

1 /a tot do/ dx F 



and W = 
defined by the condition xp > 

1 /a tot do/ dx F 




10 2 



_5_ 



(a) 



1 

Xf 



(b) 



Fig. 2. Distributions in xf in comparison with HI (a) and ZEUS (b) data. 

The parton density parametrizations are: GRV LO [ ] (with leading-order 

matrix elements); GRV HO [ ], MRS A' [---], CTEQ 3M [ ] (with 

next-to-leading- order matrix elements). 



3 We note that the present calculation does not permit a next-to-leading-order 
comparison of transverse momentum spectra. 
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The HERA proton and lepton energies are Ep = 820 GeV and E\ = 26.7 GeV, 
respectively. The cuts imposed in the experimental analyses are: (a) HI cuts: 
10" 4 <x B < 1(T 2 , Q < 10 GeV, 54.77 GeV < W < 200 GeV, E v > 14 GeV, 
157° < -& v < 172.5°; (b) ZEUS cuts: 0.04 < y < 0.85, 3.16 GeV < Q < 
12.65 GeV, 75 GeV < W < 175 GeV, E v > 10 GeV. Here x B = Q 2 /2Pq, 
y = Pq/Pl, and Ey and are the energy and polar angle of the outgoing 
lepton in the laboratory frame. The experimental data are corrected from the 
limited detector acceptance for the observed charged particle to the full phase 
space. 

Figure [| shows the comparison of experimental data with the theoretical pre- 
diction for various sets of parton distribution function^]. The employed value 
for Aqcd is the one from the parton distribution functions. Charm and bot- 
tom quarks are treated as massless flavours in the matrix element, and the 
flavour threshold in the running coupling constant is assumed to be at the 
single heavy quark masses. Except for very small and very large values of 
Xp, the agreement of theory and experiment is quite satisfactory. The next- 
to- leading-order corrections are negative for most bins in xp, and bring the 
prediction in better agreement with the data, in particular in the large- a; ^ 
range. Therefore, the QCD corrections are important to describe the data. 
For large xf, the theoretical prediction is, however, systematically larger than 
the experimental result, but still within about one standard deviation of the 
data, depending on the parton density parametrization under consideration. 
It is expected that the fragmentation function picture breaks down at small 
xf, because soft fragmentation effects, such as particle production between 
the current and remnant jets, become important. Moreover, the fragmenta- 
tion function fit in Refs. ||f|] only takes into account data for 0.1 < xp < 0.8. 
Outside of these bounds, the fragmentation functions are extrapolated. We 
also wish to note that the experimental data include all produced charged 
hadrons, whereas the theoretical prediction, due to a lack of a suitable set of 
fragmentation functions, is restricted to charged pions and kaons. 

We finish this section with some remarks concerning the scale dependence of 
the theoretical prediction. Figure ||] shows the dependence of the cross sections 
and of the charged multiplicity cr/a tot , for the cuts of the ZEUS analysis, 
as a function of the factorization scales fif and fip,, varying individually in 
the form (Q, the other scales being kept fixed at Q. A minimum value of 
2 GeV for the factorization scales is required. The parton densities are given 
by the GRV parametrizations. For the one-particle-inclusive cross section, the 
dependence on fif is much flatter in next-to- leading order than in leading 
order, whereas the dependence on fip) is slightly larger in next-to- leading order 
(although it develops an extremum). Including only the scale-compensating 



4 The parton density parametrizations are: GRV LO and GRV HO ||, MRS A' 
§ and CTEQ 3M [0]. 
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Fig. 3. Scale dependence of cross sections (a) and multiplicity (b) for the ZEUS cuts. 

The particular scale set to CQ is given by: (a) [if ( )i l^D (— ■■ ) in leading order 

and [if ( ), [ir> (— ■ — ) in next-to-leading order for the one-particle-inclusive 

cross section a, [if in leading ( ) and next-to-leading (■■ ■■ ) order for the total 

cross section a"tat>' (b) [if ( ), [id (— ■■ ) in leading order and [if ( ), [ip 

(— ■ — ) in next-to-leading order for the multiplicity. The other scales are fixed to be 
equal to Q. 



terms from the next-to-leading-order matrix element proportional to a product 
of a splitting function and a logarithm in fi D , the /independence in next-to- 
leading order is much smaller than the one in leading order. The enhanced [in- 
dependence is therefore a genuine higher-order effect. The /i/-dependence of 
the total cross section is reduced in next-to-leading order as well. Although the 
scale dependence of the absolute cross sections looks quite reasonable, the scale 
dependence of the multiplicity is larger in next-to-leading order than in leading 
order (Fig. ^]b). A similar study shows that the overall scale dependence of the 
multiplicity is significantly reduced for larger values of Q. The impact of the 
scale uncertainty on an a s -measurement will be studied in the next section. 
The dependence on the renormalization scale, not shown in the figures, is very 
small, because the running coupling constant enters the matrix elements only 
in next-to-leading order. 
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3 The Strong Coupling Constant via Scaling Violations? 



The good agreement of the XF-distributions of theory and data, as discussed 
in the previous section, is an encouraging sign for the feasibility of a measure- 
ment of the strong coupling constant via scaling violations of fragmentation 
functions. The scale evolution of fragmentation functions is governed by a 



renormalization group equation ||1 1| 



dD h/j (z,p 2 D ) 
d\n fj, 2 D 



<Xs{l4>) 

2tt 



Dh/k [ — j 



(3) 



A multiparameter fit of fragmentation functions for partons into charged 
hadrons thus permits a determination of the fundamental QCD parameter 
Aqcd (or equivalently the strong coupling constant a s (M§) at the mass of 
the Z boson), if the measurement is performed for several scales Since 
the scale evolution depends only on the logarithm of the scale, it is necessary 
to include a wide range of scales in the fit. The photon virtuality Q is the 
only scale at hand in the process under consideratioiQ, and consequently the 
analysis will require high statistics, because of the rapidly falling cross section 
with increasing Q. 

We study here three potential sources of errors: the statistical error due to 
a limited number of events, in particular at large Q; the error coming from 
the spread of the theoretical cross sections due to various parametrizations of 
parton densities, and the uncertainty due to the choice of the factorization 
scale he). The error estimate presented here can, of course, not replace an 
analysis based on a y 2 -^] it is rather intended to give a first assessment of 
the sizes of the various uncertainties. 

The procedure to obtain an estimate of the statistical error of a s is the fol- 
lowing. We fix the fragmentation functions at a scale of po = 2 GeV as the 
leading-order parametrization of Refs. @,f§0. We then evolve this input with 
two different values for Aqq D of {a} 0.1 GeV and {b} 0.2 GeV. The corre- 
sponding ^-distributions p^ a ' and are determined for these two sets of 



The evolution kernels (u,a s ) in next-to-leading order can be found in 



Ref. U; see also Ref. fll3|. 

6 The total hadronic energy W is not directly related to the hard scattering process, 
because it contains contributions from the proton remnant jet. 

7 We have to make an assumption of this kind in order to get an estimate of 
the dependence of the x^-distribution on the value of the strong coupling constant. 
Since later on we take ratios of the distributions at two different scales (we are only 
interested in the slope depending on Q), it is not really important at which scale 
we identify the distributions. 
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fragmentation functions. With the cuts of the ZEUS data analysis unmodi- 
fied, except for the cut on Q, we determine the distributions 

p {al}^ p {a2} ; p {a3} ) 

p{6i} ; p {b2} ) p {b3} for three binf g in q. {!} [3 16; 12 .6] GeV, {2} [12.6, 100] GeV 
and {3} [100, 150] GeV. The ratios A* 21 * = p^Vp {1} and A* 32 * = p^/p {2} 
for an arbitrary coupling constant a s (taken at the mass of the Z boson) are 
assumed to depend linearly on a s : 

X «J = A «J + ^^i («._„«). (4) 
a} - a} 1 



For a given luminosity, and under the assumption of a Gaussian error for 
the event numbers, this allows us to estimate the statistical error e& of a s for 
every bin k in xf- These individual errors are then combined into a total error e 
according to e = 1/ ^Y.k(^/ e l)- To obtain explicit numerical values, we use 
the CTEQ 3L parametrization [[Uj for the parton densities (for simplicity, we 



work in leading order). The integrated luminosity is assumed to be 250 pb -1 
(the HERA design luminosity, per experiment, integrated over five years). For 
the analysis based on bins {1} and {2}, we obtain a statistical error of a s (M|) 
of ±0.0007, and for the bins {2} and {3}, the statistical error is ±0.027. 

Another problem of the extraction of a s is the dependence of the theoretical 
cross sections on the parton density parametrization. Taking the ratio of one- 
particle-inclusive and total cross sections leads to a substantial cancellation, 
but a certain residual dependence remains. To estimate the size of this effect, 
we determine the spread of the results for a s (M|) depending on the next-to- 
leading-order parton densities from Refs. []S|j9|JT0[] . For the bins {1} and {2}, 
the spread is ±0.017, and for the bins {2} and {3}, the spread is ±0.005. 
Future global fits of parton densities including improved HERA data should 
reduce this systematic uncertainty. 

Finally, we discuss the dependence of a s on the choice of the factorization 
scale fi£>- To obtain an estimate, the ratios A are determined for the three 
choices Q/2, Q and 2Q of this scale. The change of cross section has for 
consequence a variation in the extracted a s (M§)-value of ±0.013 and ±0.011 
for the combinations of the bins {1}, {2} and {2}, {3}, respectively. 



It turns out that the dependence on the parton densities is unacceptably large for 
smaller values of Q due to the large spread of parametrizations at small x. Bin {1} is 
the Q-range of the original ZEUS analysis. The other two bins are chosen such that 
one of them includes data only at very large Q. We expect that bin optimization 
will lead to a smaller overall error by balancing statistical and systematic errors; 
this is, however, beyond the scope of the present paper. 



7 



4 Summary and Conclusions 



We have calculated, in next-to-leading order, the £p-distribution for the pro- 
duction of charged mesons in deeply inelastic scattering at HERA. The com- 
parison with experimental data is satisfactory. The dependence of the the- 
oretical cross section prediction on the factorization scale fif is consider- 
ably reduced, whereas the scale dependence of the charged multiplicity is 
increased. The /z^-dependence of the next-to- leading-order cross sections, al- 
though slightly larger than in leading order, is small. The cross section results 
can therefore be considered to be reliable for a quantitative test of QCD. 

As a particular example, we have estimated the statistical and systematic 
errors (due to theoretical uncertainties) for an extraction of the strong cou- 
pling constant via scaling violations of fragmentation functions. We find, for 
the HERA machine running at design luminosity for five years, Aaf at w 
±0.0007/±0.027, Aa PDF « ±0.017/±0.005 and Aaf alc « ±0.013/±0.011 for 
two specific sets of cuts. It is expected that there is still some room for im- 
provement by means of bin optimization. We have not considered experimen- 
tal systematic uncertainties. They will also contribute to the error on a s . 
The dominant systematic error in the experimental x^-distributions comes 
from the boost of the particle momenta from the laboratory system into the 
hadronic centre-of-mass system. For a measurement of a s it may therefore 
be advantageous to define a suitable observable in the laboratory system, to 
circumvent this source of systematic uncertainty. 

Compared with the present error Aa s = 0.006 of the world average, the errors 
estimated here are large. It should be kept in mind, however, that a reduc- 
tion of Aa PDF can be expected because of improved fits of parton densities, 
in particular in the small- x region. It might also be possible to reduce the 
dependence on the parton density by performing the analysis in the Breit 
frame [14]]. A reduction of Aaf ale would require a higher-order calculation or 



a physical scheme for the choice of the factorization scales. A measurement of 
a s at HERA via scaling violations of fragmentation functions is worth doing 
because it is an independent quantitative test of QCD and, more important, 
because it complements other methods for an a s -determination in deeply in- 
elastic scattering, such as the measurement via (2+l)-jet rates and via scaling 
violations of structure functions. 
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